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Intensity mapping of the neutral hydrogen brightness temperature promises to provide a three- 
dimensional view of the universe on very large scales. Nonlinear effects are typically thought to 
alter only the small-scale power, but we show how they may bias the extraction of cosmological 
information contained in the power spectrum on ultra-large scales. For linear perturbations to 
remain valid on large scales, we need to renormalize perturbations at higher order. In the case 
of intensity mapping, the second-order contribution to clustering from weak leasing dominates the 
nonlinear contribution at high redshift. Renormalization modifies the mean brightness temperature 
and therefore the evolution bias. It also introduces a term that mimics white noise. These effects 
may influence forecasting analysis on ultra-large scales. 


I. INTRODUCTION 


Mapping of the integrated 21cm emission from neutral hydrogen (HI) within unresolved sources will allow us to 
cover very large volumes of the observable universe with upcoming radio telescopes. This promises to deliver excellent 
precision in measuring the standard cosmological parameters [T], as well as new opportunities to probe signals on 
ultra-large scales, such as primordial non-Gaussianity and horizon-scale relativistic effects [H]. 

The ability to extract precise cosmological information from an experiment is dependent on accurate understanding 
of the theoretical model. On scales k < k^i = 0.2/i(l -|- Mpc“^ (where = 0.96 is the primordial spectral 

index), gravitational evolution of large-scale structure can be analyzed via cosmological perturbations [TU]. For scales 
k ^ fcni, it is commonly assumed that one need not go beyond linear order. This is valid for cold dark matter 
(CDM) perturbations, which obey precise energy and momentum conservation, but not necessarily for tracers of the 
underlying CDM distribution. 

It was shown in m that mode coupling in the biased density field contributes power on large scales which becomes 
important on Hubble scales. This potential threat to the validity of linear perturbation theory on large scales can be 
mitigated by renormalization [HHH]. Renormalization leads to a change in the mean (zero-order) density and in the 
bias parameters for the tracer. In addition, it introduces a noise-like term [I1II3II1]. 

We employ this renormalization technique to study the HI brightness temperature after reionization. Up to second 
order in perturbation theory, the mean HI brightness temperature is modified by a term that depends on redshift and 
a coarse-graining scale. At z < 1.5, this term is proportional to a nonlinear bias parameter, while for z > 1.5, it is 
dominated by a nonlinear weak leasing term. We point out that the renormalization of mean brightness temperature 
leads to a renormalization of the evolution bias parameter, which modifies the linear power spectrum on horizon 
scales. The main implication is that cosmological constraints on horizon scales need to incorporate nonlinear effects 
that flow from renormalization. 

Our fiducial model is determined by the Planck 2015 best-fit values [H!; in particular, h = iJo/100 = 0.673, 
^mo = 1 ~ ^Ao = 0.315. A b usy reader can skip section |H] which covers the perturbation theory (with details given 
in Appendix A). Section HI covers the nonlinear (Gaussian) bias (with details for HI in Appendix B), and the basic 
ideas of renormalization are introduced and then implemented in redshift space. We discuss our results in section IV 
and conclude in section El 


II. FLUCTUATIONS IN HI BRIGHTNESS TEMPERATURE 


The brightness temperature observed (in redshift space) is given by [ 2 ] 

T°‘’"(z,e) = ^|-^n(z,e)J(z,e), 
4 /csUai 


( 1 ) 


where z is the redshift of the source, e is the unit direction of the source, E 21 is the proper energy of the emitted 

3-15s-1 

of the map between real and redshift space [see (A18)]. 


photons, and Aio = 2.869 x 10 ^°s ^ is the emission rate. The number density of HI atoms is n and J is the Jacobian 
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The observed temperature is expanded up to second order: 

e) = f{z) [l + (z, e) + ^ (z, e) 

where we use At to denote the quantity observed in redshift space. The mean brightness temperature is 


T(z) = 


Stt^ hz'AiQ n{z)a{zY 


4 kBE 2 i H{z) 


+ mK. 


0.003 


H{z) 


( 2 ) 


(3) 


Here fini is the comoving HI mass density in units of the current critical density. 

The observed fractional temperature perturbation At is related to the fractional number overdensity Sn in real 
space as follows (see Appendix A for details): 
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(5) 


Here fc“ = da;“/dA is the photon 4-momentum, Aa;“ is the perturbed physical position, oj is the metric vector 
perturbation and || and T denote projection along and transverse to e. The evolution bias is defined by the mean 
number density: 


be = -y{z) 


dln[n(l + 2 ;) 


dln(l + z) 

From now on, for convenience we will often drop the superscript (1) for first-order perturbations. 
• Background: We calculate Hhi from a simple halo model using 

1 pM 


( 6 ) 


^hi(.2:) — 


(1 -7 zY Pco 


with Phi ( 2 ) = 


dMn^(M)MHi(M), 


(7) 


where pco = 377g /(SttG) and Uh is the halo mass function. See Appendix B for details. 


First-order perturbation: In Poisson gauge [see (All], the full first-order perturbation is given by (A30|: 

77' 


A^’(z,e) = (5);’ - -7 ^6, - 277 


77^ 

y 


" 77 77 


he-m- 


77 


$ 


( 8 ) 


1 


77 


77' 


+ - 6e- 277 - — / dAC^'T^-'), 


77 


where v is the velocity potential. We h ave dr opped the subscript s for the source and set the observer-dependent 
monopole and dipole part to zero [see (A30)]. Equation (|^ is in agreement with (17) of [2]. 

The second term is the linear Kaiser redshift space distortion (RSD) term, followed by the Doppler term, the 
rate of change of the curvature potential, the gravitational redshift term and finally the integrated Sachs-Wolfe 
term. The third to the last terms are due to relativistic lightcone effects, and are suppressed on subhorizon 
scales [2]. Note that there is no weak lensing contribution to the HI temperature at hrst order. 

A further horizon-scale term arises when we introduce the scale-independent bias. The relation Sn = bSm, where 
b = b{z), is only physical (and thus gauge-invariant) if expressed in the matter rest-frame, i.e. in comoving- 
synchronous gauge [2]: 


6(^1) = + (377 - be)v. 


(9) 
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• Second-order perturbation: The second-order perturbations of the galaxy number counts that are observed 
on the lightcone in redshift space, were recently computed by [18H22) . The extremely lengthy expressions are 
however not directly applicable to the case of intensity mapping which we are treating. In Appendix A we 
compute A^’(z, e), using the approach developed in [531 [H]. The full expression is unmanageable, and we make 
a consistent approximation that generalizes the Kaiser approximation [25] to second order: 
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nonlinear lensing 


nonlinear time perturbation 


The nonlinear RSD terms recover the result of mills] for the galaxy number density. The nonlinear lensing term 
is i.e. the contraction of the weak lensing deflection angle with the transverse gradient of the 

first-order (density + RSD) term. This is the dominant nonlinear correction to [111I26j . The time perturbation 
term is subdominant in the power spectrum - but it may be more important in the bispectrum. There are many 
further terms that we have neglected because they are strongly suppressed on sub-Hubble scales. 


Relationship with galaxy number density 


We derive a covariant mapping between the perturbation of observed brightness temperature and the perturbation 
of observed galaxy number density. If we count the observed number of HI atoms per solid angle and per redshift 
instead of galaxies, we find that the observed brightness temperature is given by (5] 


_ h^Aio n°^^(z,e) 
4 Da{zY 


( 11 ) 


where Da is the angular diameter distance. Expanding in perturbations we find 
JDa 
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.(13) 


Here A„ is the fractional HI galaxy number overdensity observed in redshift space, before magnification bias is applied. 
At linear order, there is an accidental symmetry whereby A^^ is obtained from the magnification bias corrected A„, 
by simply setting the slope of the luminosity function to s = 2/5 [5]. Equation (12) is in exact agreement with (23) 
of |5] and (131 is its second-order extension. 


III. NONLINEAR PERTURBATIONS AND RENORMALIZATION 


In order to illustrate the main idea, we start by ignoring the corrections due to observation in redshift space and 
by focusing on sub-Hubble scales, where the difference between comoving-synchronous and Newtonian gauges in © 
may be ignored. Fluctuations in number density should average to zero over large enough volumes, i.e. 


^n(^,x) 


n(z,x) - {n){z) 
{n){z) 


((5„(z,x)) = 0. 


(14) 


Within standard cosmological perturbation theory, we normally assume that {n){z) = n{z) = h^'^^'^( 2 :), i.e we assume 
that fi is determined by the background FLRW spacetime. However this does not take account of the complexities of 
relating a tracer to the underlying CDM distribution. 

A bias model is needed to relate Sn to the CDM overdensity Sm- We assume that the primordial perturbations are 
Gaussian, and use an Eulerian local bias model, applied up to second order: 

1 2 

^„(z,x) = bi{z)6m{z^^) + 2^2(2) [5m(2,x)] . 


(15) 
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As a consequence, the average in (15) no longer vanishes: 


{5n){z) = b 2 {z)al{z), al(z) = f 

fcmin ^ / 


(16) 


Here tT| is the variance of Sm and ks is the small-scale cutoff. 

A mathematically consistent way of ensu ring zero average is via renormalization [12111811151116| . This involves 
subtracting 620'|/2 from both sides of ( [l^ . The right-hand side will then have zero average, by (16), while the 
left-hand side will define a renormalized nuctuation = 6n — b 2 crg/ 2 . 

Renormalization of Sn must leave invariant the physical number density n. This is achieved by renormalization of 
the background number density, n, which in turn induces a renormalization of the bias parameters: 




n 


— n{l + b 2 ag/ 2 ), b^— 


(1-k 62 cr|/ 2 ) ■ 


Putting this together, we have 


= 


= b^Sm + 7.^2 “ Cr|] 


(0 = 0 . 


(17) 


(18) 


Note that the fundamental physical quantity (i.e n) is independent of the cutoff scale chosen to regulate the integral 


in tT|, i.e., dn/dks = 0 [IB] . 


An important consequence of renormalization was not mentioned in [13 nans], since these papers did not include 
the horizon-scale general relativistic terms that we have in (|^. This consequence is that the evolution bias bg, which 
appears in the coefficients of most of the general relativistic terms, will be renormalized as a result of renormalization 
of the mean number density. Using equations ([^ and (0, be is modified by n-renormalization as 

dln(l -I- 620-1/2) 


bR = be-n{z)- 


(19) 


dln(l-hz) 

As a result, the linear power spectrum of the observed fractional temperature fluctuations will be modified on horizon 
scales by renormalization. 


The power spectrum may be computed by expanding (18) in Fourier space (we suppress the z-dependence): 


6^(k) = 6f<5„(k) + 


1 




,(20) 


^^:^^^;^^m(ki)6„(k2)(27r)36^(k-ki - k2) [6f-t 6f F2(ki, k2)] - 6f(5^(k) 

where Pm is the linear matter power spectrum and F 2 is the standard nonlinear density kernel, which we approximate 
using a matter-dominated model: 


F2(ki,k2)=^ + ^ 
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k2 
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ki ■ k2 
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( 21 ) 


Computing (i5,(^(k)(5^(k')), we find that the contributions of the a-| term exactly cancel out and the power spectrum 
of the renormalized fractional number overdensity is 

d^ki r, 


Psn{k) = {b^) Pm{k) + 


1 2 


(27r)3 I 


6" + 6fF2(ki, k - ki) P^(|k - ki|)P™(fci). 


( 22 ) 


(23) 


The b 2 term leads to constant power on very large scales [ni[I3l27]: 

k^O ^ PsR{k)^-{b2) J 

where we used F 2 (k, —k) = 0. It is possible to renormalize the power spectrum by subtracting this constant power 
from both sides of (22). In this approach, the constant-power contribution is treated as part of the noise budget 

[nHI3: 

P«(fc) = (6f)'P™(fc) + i I ||^|[6f + 6fP2(ki,|k-ki|)]V^(|k-ki|)-(6f)V„(A:i)}p„(fci), (24) 

Psnik) = p,l{k) + iVeff, iVeff = ^ (6^)' J ^P^{k,). (25) 

Here we follow a different approach, consistent with im nHj. We see the constant-power term as containing 
cosmological information rather than noise, i.e. we treat it as a nonlinear correction to the power spectrum. In this 
view, we use the unrenormalized power spectrum (22). The two approaches differ only in the interpretation of the 
correction to the linear power spectrum on ultra-large scales. 
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Renormalization in redshift space for the brightness temperature 


We now extend the established renormalization procedure from real space to the observed redshift space, in the 
case of the observed fractional HI brightness temperature perturbation: 


Ariz^e) = 


' (r°bs)(2) “ 


(26) 


The sky average of (Ap) is zero by definition, but if we assume that gg^ nonzero value. In 

order to ensure that (A^) = 0, we have to renormalize the background temperature. 

In observed redshift space, the average is more complicated because of the nonlinear leasing contribution in (101. 
After performing the sky integral, we find that 
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(28) 


is a transfer function for the weak leasing convergence. The sky average of the nonlinear leasing term is the only 
redshift-space contribution to the monopole. 


We renormalize the background temperature and its fluctuations by adding and subtracting (At) to (26): 
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where we suppressed the z-dependence. Here the renormalized first- and second-order parts are 
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(29) 

(30) 

(31) 

(32) 


In order to define the HI bias physically at linear order, we used the definition (|^ and we applied © in comoving- 
synchronous gauge. For convenience we have omitted the cs superscript - i.e., 6^ in the above and following equations 
is understood to be S^. Note that we do not need to extend the comoving-synchronous definition to second order, 
because the difference contains only terms that we are neglecting at second order. 

The A and B terms in (31) are the coefficients of the first-order potential and Doppler terms in (1^: 
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(33) 

(34) 


The renormalized background temperature is equivalent to a renormalized background number density by ([^ , so that 
the evolution bias is renormalized: 


b^ = b,-n{z) 


din (1 -I- {At)l) 

dln(l -I- z) 


(35) 


We have introduced a smoothing scale L (in units of h ^Mpc). For a top-hat window function on a scale L, (27) 
becomes: 
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(36) 


We show the fractional difference, Ab^/bf, = {bf — 6e)/&e, in Fig. for various values of L. For L = lO/i ^ Mpc, 
the difference is ~ 1% and for L ^ 10h“^Mpc, we recover the FLRW limit as expected. The difference grows as 
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L is decreased below lQh~^ Mpc. For L < 2h~^ Mpc, we need to go beyond second-order perturbation theory. The 
crossover point at z ~ 1.5 arises due to competititon between the nonlinear HI bias and the nonlinear leasing effect. 
The nonlinear bias in (35) is negative and the leasing contribution is positive, and they cancel at z ~ 1.5. The 
fractional difference is dominated by the nonlinear bias for z < 1.5 and for z > 1.5 by the nonlinear leasing effect. 



FIG. 1. Fractional difference in the evolution bias due to nonlinear effects for different values of the smoothing scale L (in units 
of Mpc). 


The kernels in (32) are as follows. G 2 (ki, k 2 ) is the kernel for peculiar velocity at second order, which we approximate 
by its form in a matter-dominated universe: 


G2(ki,k2) 


7 2 \k2 kj fcife 7 'v fclfc2 J 


(37) 


The /C kernels encode the nonlinear effects from mode coupling: JC^ - RSD (velocity-velocity and velocity-density); 
/Cl ” leasing (lensing-density and lensing-Kaiser); /Ct - time perturbation (velocity-density and velocity- velocity). 
These are given by: 
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The power spectrum of observed HI temperature fluctuations in a distant observer approximation then follows from 
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(31) and (32): 
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Once again, this power spectrum has a nonzero constant-power limit as fc —> 0: 
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The monopole of the power spectrum is 
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where Pi is the second order part of (42). 


IV. RESULTS AND DISCUSSION 


There are two key nonlinear effects on the ultra-large scale power spectrum that we have identified - modification 
of the evolution bias and constant power that mimics white noise on horizon scales. These could potentially affect 
the standard analysis on horizon scales which assumes that linear theory is valid. 

Here we compute the nonlinear modifications to the power spectrum Pt of the renormalized HI temperature 
fluctuations, as given in (42) and (44). For this computation, we compute the HI bias parameters using a standard 
prescription in Appendix B. We set the limits of the convolution integral in (42) following [21]: (fc“ 
(10-4,104)/i-iMpc. 


(^min 7.,max 


) = 




FIG. 2. Power spectrum Pt of the renormalized observed HI brightness temperature, at z = 0.5. Left: along the line of sight; 
Right: transverse to the line of sight. The solid black curve shows the total power and the red curve is the linear part of 
(42 I. The remaining curves show the nonlinear effects from: F 2 - second-order density perturbations, G 2 - RSD from second- 
order velocity, 62 - HI bias, /Cr - RSD from mode coupling, /Cl - lensing effect from mode coupling, /Ct - time perturbation 
from mode coupling. 


In practice, the choice of the smoothing scale \/ks will be connected to the resolution of the experiment. This 
scale cannot however be smaller than the scale I/fcni if we wish to remain in the regime of validity of the perturbative 
approach. In our calculations we set ks ^ k^\ L~^. We note that although the results can be sensitive to the exact 

value of fcni, this is not an issue as long as I/Zcs > I/fcni- The calculations will then only depend on ks but this is 


































just a statement that the power spectrum needs to be smoothed on the resolution scale of the experiment in order to 
compare to observations. 

We do not use a window function to regulate the convolution integral in (42). Careful analysis using a smoothed 
density field and angular power spectrum will be presented elsewhere. 




10 "^ 10 "^ 10 "^ 10-1 10 " 
i:[hMpc-i] 


FIG. 3. Power spectrum Pt of the renormalized observed HI brightness temperature at various redshifts. Left panels: along 
the line of sight; right panels: monopole. Upper panels: Solid lines are the linear power spectrum, dashed lines the nonlinear 
power. Lower panels: The fractional difference due to nonlinear effects. Vertical lines indicate the comoving Hubble scale. 


Figure illustrates the amount of power contributed by the terms in (42), for fj, = 1,0. The term responsible for 
non-vanishing of the average of contributes constant power on large scales. At linear order, the horizon-scale 
GR effects in become important near the horizon (fc ^ H) and dominate on super-Hubble scales {k <%). (We 
have not computed the horizon-scale GR effects at second order, since this in itself is a major task which is left for 
further work.) Along the line of sight [pL = 1), the nonlinear RSD from mode coupling (via /Cr) clearly dominates 

the nonlinear contribution at fc > 0.05. For k > 0.1, it dominates over even the linear power spectrum P^\ At 
BAG scales, its contribution is ~ 10%. As expected /Cr vanishes for fj, = 0. 

We quantify in Fig. the total second-order effect on the linear power spectrum, focusing on large scales. The 
fractional difference is defined as APp/P^^ = [Pp — Pp^]/Pp''- The contribution of the dominant constant power 
term peaks around the horizon scale k = 'H(z). The constant power from the lensing term /Cl dominates over the 
constant power from nonlinear bias at about z ~ 1 . 


clT 

<1 



FIG. 4. Effect of renormalized evolution bias on the monopole of the power spectrum on large scales. (The line of sight power 
spectrum is identical.) The smoothing scale for the evolution bias is set at L = 5(1 -I- Mpc. 
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We also quantify the error on the power spectrum associated with using the bare evolution bias ([^, instead of the 
renormalized evolution bias (35). To this end, we define the fractional difference as APt/Pt = {Pt — Pt) /Pt, where 
Pt is the power spectrum with bare evolution bias, given by (j^. Figure]^ shows that the difference could be more 
than 20 % on horizon scales. 

Furthermore, what about third-order terms that we have neglected? These terms can contribute to the one-loop 
correction to the tree-level power spectrum through cross-correlation with first-order terms. Here we argue that 
the effect of third-order terms will be sub-dominant on the large scales that we are interested in. The brightness 
temperature up to third order is 


T°'°\z,e)=T{z) 


1 +A: 


T + o^T (^.e) 


A^>(z,e) 


(45) 


Assuming Gaussian initial conditions, the main contribution to the mean brightness temperature after renormalization 
will come from the third-order bias term. For a local Eulerian bias model, = bi6m + b2S‘^/‘2 + b3S^/3\, the ensemble 
average of the third-order term is given by m- 


63(e) = 3&3((0'0 = 


(46) 

Using a Sheth-Torman mass function [32l |33], we can show that —1 < 63 < 1 for z < 3.5, as seen in Fig. (see 
appendix [^. This implies that for cr| < 1, we have (cr|)^ < ^s- Therefore (A^') < (A^’), i.e the contribution from 
third-order fluctuations to the mean brightness temperature is sub-dominant. 

The mode coupling term (A^’(z, k)A^'(z, k)), will contribute to the power spectrum at one-loop as follows: 


3 ( 13 ) 


(fc,M) = ^ [bi+ +a'^ 
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Pmik) 


d^fci 

(27r)3 


P„(fci)/C<^>(^,ki,-ki,/r), 


(47) 


where is the kernel for all third-order terms. Although (47) is a one-loop contribution to the tree-level power 


spectrum, it does not appear as a convolution of two power spectra, i.e. it does not produce a coupling of power on 
small scales to ultra-large scales that can modify our results on those large scales. Also if we set ki = kr, and take 
the limit fc —>■ 0, the pre-factor in (47) goes to zero faster than the rate at which the r-integral goes to larger values. 

Finally, it is important to point out that the expression for HI brightness temperature given in equation ([^ is a 
special limit of the HI brightness temperature during the Epoch of Reionization (EoR): 


ST = f Xhi (1 + 6 „) ( 1 + ^ 


‘"t) 


(48) 


where Xm is neutral fraction of hydrogen, (5a;Hi is the perturbation in neutral fraction of hydrogen, Tcmb is the 
temperature of the CMB, SJ is the perturbed Jacobian and Tg is the spin temperature. Equation (48) reduces to 
equation ([^ in the limit Tqmb Tg. Despite these differences, the modulation of the HI power spectrum on ultra- 
large scales by nonlinear effects we discussed at low-z will also be applicable during the EoR because the dominant 
nonlinear terms come from the perturbed Jacobian (SJ) and not from Tg and Tcmb which we neglected. 


V. CONCLUSION 

We derived for the Hrst time the second-order perturbation of the HI brightness temperature in the post-reionization 
universe, focusing on the dominant nonlinear redshift space distortion terms and the nonlinear lensing term. We find 
that: 

• The physical mean density of HI or the mean brightness temperature is modified by nonlinear terms that depend 
on redshift and the cutoff scale. By smoothing over structures on scales where perturbation theory may not be 
trusted, we showed that this could lead to more than 20 % change to the power spectrum on large scales. 

• The constant power on large scales associated with nonlinear evolution of biased tracers can significantly affect 
the power spectrum on horizon scales. The boost in amplitude could affect the signal to noise ratio for possible 
detection of general relativistic effects and primordial non-Gaussianity from intensity mapping. 

However, the power spectrum is not observable. In order to determine whether these effects do lead to observable 
consequences, we need to compute the observed angular power spectrum C^(z,z'). This is a subject for further 
investigation. 
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Appendix A: Brightness temperatnre in pertnrbation theory 


Here we provide details for the perturbation calculation up to second order, using the approach developed in |231l24j . 
We consider a perturbed FLRW spacetime in the Poisson gauge assuming a flat background metric 

ds^ = a^ds^ = [—(1 + 24) + + 2wida;*d?7 + ((1 — 2'1> — + hij) dxMx-^] . (Al) 

Here the proper time is related to conformal time by dt = adrj. $ is the first-order Newtonian gravitational potential, 
'k is scalar curvature perturbation, and uJi = and hij = are the second-order vector and tensor contributions, 
their divergences vanish. Since the null structure is unchanged under conformal transformation, the metric gab 
associated with the line element ds^, i.e physical spacetime maps the null geodesic of the physical spacetime gab 
to a null geodesic on a perturbed Minkowski space time gab with the affine parameter associated with each metric, 
transforming as dA —>■ dA = a“^dA. The photon 4-vector transforms as ^ ka = ka- For the matter 

4-velocity, we have u°‘ = a~^u‘^ Ua = aua ■ Hence, the photon energy transforms as E = —Ubk^ = —aubk^ = aE. 
For any tensor S we expand it up to second order as 

S = S + 5^^^S+]^5^^'>S. (A2) 

We expand the 4-velocity, of a matter field using ( |A1[ ) up to second order 

yO = 1 - $ -h -h ]^d,vd^v 

-k 

The perturbed photon tangent vector is calculated from the geodesic k‘^\7ak^ = 0 and the redshift is given by 

l + 2 = l + (A5) 

The background conformal time or the conformal metric affine parameter is mapped to the conformal time in redshift 
space according to 

A = ^-k<5A-ki(5<">A (A 6 ) 


(A3) 

(A4) 


Hence the physical scale factor is related to the background scale factor through Taylor series expansion 


a(A) = a^Vs) 


1 + HsSX + ^ {Hs <5'">A + (H( + Hi) (<5A)2) 


(A7) 


where we have performed the Taylor series expansion around the position of the physical position of the source. For 
consistency, we could define a homogeneous redshift in terms of the physical scale factor as 


(1 -k Zs) — 


aj’^o) 

a(^s) 


We have also made the following replacements in (A71 


1 da 
a dA 


= n - — 

adA2 


= (K + nl) ., 


(A8) 


(A9) 


Inverting (A51 and using (A7) leads to 
1 a{i^s) 


(1 -k Zs) a{vo) 


1 -k {Hs5\ - + 2 (77s5'"’A - -k 2{5z)'^ - 2HsSz5X + (U's + Hi) (JA)^) 


(AlO) 
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Thus from (AlO), the background scale factor may now be re-defined in terms of the physical redshift using (A 8 | and 
requiring that terms vanish at every order leads to 


At 






(All) 

(A12) 


the affine parameter associated with the physical metric according to {drj —dA) 


For further details on this see [24]. In order to expand the Jacobian in we need to relate the conformal time to 

(A13) 


dA 


dA . 


1 


l + dk° + - (W|| 


and using (A5| we link it to the physical redshift 


dA , /v \ 2 

= (1 + ^s) 

dA 


1 


1 -b {Sk° - 2Sz) + - (W|| -b + 2Sz (3<5z - 2Sk°)) 


Differentiating (A51 wrt A gives 

dz 'Hiv) 


1 , , 1 / dJz 1 d5< 

l + 5z+ -5^^z - —— -r— -b - 

2 '^( 77 ) \ dA 2 dA 


dA 0 ( 77 ) [ 

It is convenient now to re-express background ^( 77 ) and 0 ( 77 ) in terms of the observed redshift 


F - 

dA a{v) 


1 -b 


1 dSz 

f n'd 


1 

dz — — -r-- 

[ns-^ 

1 JA 


ns dA 

V nsj 

/ 

2 


-b2(5A 


/ dSz 
\~dX~ 


n: 




SziUs- + {5Xf + ^ 






1 dS^^^z 
W, dA 


ns 


Using (A14|, we substitute for dA and after some simplification we find 


^ _ njiy) 

dA a(7/)3 


/ n’'^ 

\ 1 ’ 

[ns- 

1 <5A 1 -b - 

V ns, 

1 J 2[ 


- - 


K 

ns 


1 d5z 
^ dA 


,5<">A + 2 ( {5zf - ^Sk°^ ) + 26z 

ns dx j 


2 d(5z 




fdSz 


1 d(5("U 
Ws dA 


— Sk I -b 2 (ns — 


nfs 

ns 


n' 


ns 


- 26 k° - XT pA + {sxf + -y 


'17 7 
IL a 


ns 


The Jacobian in 0 is then obtained by simply inverting ( |A17[ ) 

n' 


j = 


dA 

dl 


a(77)3 


n{v 

+ (ns- 


1 -b (^(5z-J/fc°-b 

n' 


ns 


<5<">A-b JA 


(JA)^ [nl-n's + 2 


1 dSz 
Ws^ 
dSz 
liX 

n'A n 


ns-^]6X]+- - 6<-^>k° - cu„ -b 

1 - ] + 2 Sz{ns- - 26 k° (ns - 


ns 


ns 




- 2Sk° Sz -b 


ns 


1 d(5<"U 
tTs dA 

n/s ^ 
ns 


1 dJzA 


ns dA / 


+ 2 {5ky + 


1 d6z 

Ws^ 


(A14) 


(A15) 


- ( - ^ ) <5<^>A (A16) 


(A17) 


(A18) 


For the HI number density we follow the standard prescription given in (A2)to perturb it up to second order 

77(77, = "^(77) 


1 




The background number density is then mapped to its equivalent in redshift space 

n{vi) = n{vs) 1 + 4^A-b ^ f 4 <5'^'A-b F((5A)^ 

n 2 \ n n 


(A19) 


(A20) 
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The rate of change of background number density along the line of sight define parameters that quantify our insufficient 
knowledge of the background spacetime, i.e the evolution bias parameters: 


be 

dbe 

Jx 


dln(a^n(7?)) a' n' 

dr] a h ’ 

- [be - m) {m - be) + m' + 


(A21) 

(A22) 


Putting (A17) and (A19) into we arrive at equations Q and ([5]). 

Using the photon geodesic k°'Vak^ = 0, we calculate the weak lensing terms from the perturbed photon 4-vector, 
which at first order is given by 


5k^{Xs) 

dk\\{Xs) 

<5fcl(A«) 






- $o) -f - To) - y ($' + T')dA 
T)dA, 


(A23) 

(A24) 

(A25) 


where we have set the perturbation of the geodesic at the observer to zero, Sk°'{Xo) = 0. The radial component is 
obtained by projecting on e*, while the transverse component comes from projecting with the screen space metric, 
Nij = gij + UiUj — eiCj. The perturbed position of the photon at first order is obtained by solving the following 
equation 


d6x^ 

~dX 


= 5k^ - 5k<^dex^. 


Using equations (A23|-(A25) in (A26) we find 
Sz 


6 X = 


1 


ns ns 


(9„ Vs - a, Ve) - (T. - To) + / (T' + T')dA 


Axy — — (To -|- To) (Ao — As) + (T -|- T) dA — —— 

J\. 


(a Vs - a,, Vo) - (Ts - To) -f / (T' T')dA 


Ax\=- (As - A)(VYT + VlT)dA = / (y - x«)(VlT + VlT)dx. 
J\o Jo 

Then the first order brightness temperature fluctuation Q is given by 


ns 


ns 


n' 


A^’(z,e) =5n- yra.Us + ^ [be- ‘AHs - yr a ■(;« + “ yr - SHs - yy T 


ns 


ns 


ns 


n's 


ns 


(A26) 

(A27) 
, (A28) 
(A29) 

(A30) 


1 

'ns 


be - ms - 


n' 


(T' + T')dA + A^>(z,Xo) 


nsy JXs 

where A^’( 2 ;,Xo) is a collection of all the terms measured at observer 


1 

ns 


Xs 

n. 


A^’(z,Xo) = [be - ns - ] d Vo + [be - ms - ) To 


1 

Hs 


T^s 

n, 


(A31) 


Recently, the second-order large-scale perturbations to the observed galaxy number counts have been computed by 
[THH52] . These results do not cover the case of HI intensity mapping. Substituting for all the terms in the second-order 
fluctuation ([^ leads to an unmanageably long expression. We take a consistent extension of the Kaiser approximation 
pS] from first to second order, so as to reduce the complexity of the calculation. The full calculation will be presented 
elsewhere. Obtaining a consistent Kaiser limit up to second order requires that we set 

6k°^0, Sk „— ^0, — ^0, 5z — ^0, (A32) 


These terms contain a Doppler term and gravitational potential terms, and products of these two terms. We do not 
set these terms to zero in the post-Born correction terms because post-Born correction terms already contain higher 
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number of derivatives than other terms in the equation. Implementing the approximation of (A321 in equations Q 
and ([^, we find 


e) « + 




1 A5z 

1 

n dA 


1 d5z 


2,5„ 


1 d<5z\ 

) 


6 z 


A^j +Ax„a„A^> + AxlV^,A^^^ 


(A33) 


. (A34) 


The last set of terms in (A34) clearly show how the effect of weak leasing appears in an HI intensity map. Its 
appearance is solely due to the post-Born correction that must be included if one goes beyond linear order |30j . The 
derivative of the redshift with respect to A contains the redshift space distortion terms and they are given by 


d6z n / -.2 r. ^ 

=5||^ -5|T + 5||® + ^ « 


-dy^\ 


(A35) 


where we have implemented the derivatives of the conditions listed in (A32). 

Considering only the dominant terms in each case and putting these terms together, we obtain finally the key result 
of (pi). 


Appendix B: HI bias from halo bias 


We calculate the bias parameters from a simple Sheth-Torman mass function [32l [33] : 


bi 

^2 

bs 


1 




2 p 


21 


Scil + {qi')P) / M 

51 (1 -f (pP) 


8., y d{p^+ vpq) -{qv-l){l + {qvY)-2p ^ P, ^2 n 

(^1 - 1) + ( —- .2.1 , ^ ((P - - 1) 


M 


236 


13 




(3 -f ivq + 3v^q^ - v^q^) 


<5? 


(8p^-I-12p^ (1-I-..g)-l-p — 2)) ^ ^ (l-|-2.2g — y + vpcYj 

H i~n / ^ ^ I ^ ~f yr, T 6 24- 


-4 


(53 (1 -I-1 -I- (i2g)P) 

(1 - ^9) , o V 


51 


51 (1 -k {vq)v) 


<53 


51 (1 -k (gi/)P) 


M 


(Bl) 

(B2) 

(B3) 


where g = 0.707 and p = 0.3 are obtained from a fit to numerical simulations. The peak height is related to the 
variance a{M) of the density field smoothed on mass scale M hy v = {5c/ uq)^ . The averaging notation is defined by 


A(z,x) 


{Xh{zy)yi = 


jy dM [Xh{z,Ti, M)Mtii{M)nh{z,yi, M)] 
jy dM [M„,{M)nh{z,^,M)], ' 


(B4) 


where M_ and M_|_ are the lower and upper limits of masses, which are related to the limits of circular velocity of 
galaxies that could house HI. These are obtained from the circular velocity constraint 


= 30\/l + z 


/ M 


VlOiOM© 


1/3 


kms 


(B5) 


where we assumed that only halos with circular velocities between 30 — 200 kms”^ are able to host HI. This range of 
circular velocity is motivated by observation [1]. 

We adopt a fitting function based on simulations for the HI mass function [4] 

Mhi(M) = CM°®, (B 6) 


where the normalization factor is chosen to match the measurement of Bhi at z = 0.8 [34]. In general the local number 
of halos with mass M is given by 


nh{M) 


= ^/(^) 


p din .2 
d in M ’ 


(B7) 
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FIG. 5. The HI bias parameters used in the paper. 


where the peak height u is related to the non-gaussian dark matter variance, anc^ ^ = (^c/cJ'nG)^ and 6c = 1.686 is 
the threshold of linear density contrast. We show in Fig. [^the shape of each of the bias parameters given in equations 
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